Flow at AGS Energies: A Barometer for High Density Effects? 
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Abstract 

Preliminary data on transverse energy 'flow' and event asymmetries reported by the 
E877(814) collaborations are compared to ARC model calculations for Au+Au at full 
AGS beam energy. ARC triple differential cross-sections for protons and pions are 
presented. Proton flow is produced in ARC, with the maximum (-Pr) ~ 120 MeV/c. 
For central events {Px) for the pions is near zero, consistent with experiment. Pion 
{Px) opposite to the nucleons' is evident in a peripheral sample, however, suggesting 
that this pion 'anti-flow' involves absorption by 'spectator' matter. 'Squeeze-out' of 
protons in central events at mid-rapidity is suggested by the ARC distributions. 
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It is presently of great interest to know how model calculations compare with data on 
flow at the AGS, especially since the excitation function of flow may well reveal new physics, 
if such is occuring at high baryon density and temperature ARC model calculations 

have previously been compared to Bevalac data on 'Sideways' flow. Other aspects of flow, 
including the so-called longitudinal and radial flows, have also been extensively studied at 
Bevalac energies. Our study of Bevalac data was done at beam momenta in the range p = 
1-2 GeV/c, for the heavy system Au+Au, energies an order of magnitude less than those for 
which the cascade was designed. Yet, very good agreement with observed flow was obtained, 
using a physically plausible treatment of two body scatterings. 

Here we are concerned with recent measurements carried out at the AGS [0], by the 
E877(814) collaboration. For details of the experimental apparatus and other information 
we refer the reader to . For our purposes it suffices to say that the instrument includes two 
highly segmented calorimeters, capable of measuring the energy deposition at various polar 
and azimuthal angles, as well as a forward spectrometer capable of particle identification. 
We will discuss calorimeter data on the azimuthal asymmetry of transverse energy as well 
as preliminary indications on the directed transverse momentum < > measured in the 
spectrometer, using the reaction plane determined from the calorimeters. 

The two calorimeters are TCAL and PCAL (target and participant calorimeters), with 
pseudorapidity coverages of —0.5 < t] < 0.8 and 0.83 < r] < 4.5, respectively. Only PCAL 
is used in constructing the transverse energy distributions. The dispersion in the event 
plane determined from the PCAL is on the order of 60° (FWHM) from experiment 0, and 
theoretical simulation using ARC events yields a similar result. Therefore relatively large 
corrections to the raw experimental numbers are needed when extracting the triple differen- 
tial cross-sections. Detector acceptances produce additional uncertainty. Nevertheless, the 
preliminary results clearly demonstrate the existence of proton fiow at AGS, and give some 
hint as to the quantitative amount. For pions one observes, interestingly, no fiow in a central 
sample 

Centrality cuts are placed on data and on ARC events using the total PCAL Et- Central 

2 



events produce more transverse energy than peripheral. For ARC we calculated PCAL Et 
using just the angular acceptance, not accounting for any efficiency corrections. The best 
procedure, of course, would be to pass ARC events through the experimental simulation 
of the PCAL response. However, the agreement in the Et spectrum calculated directly is 
very close, certainly at better than 10% level for all but the most peripheral events. For the 
purpose of making centrality cuts on the flow measurements we have assumed that E877 Et 
and (raw) ARC Et are the same. 

The azimuthal transverse energy distribution from E877 has been analysed 0, using an 
expansion in spherical harmonics. Moments Vn of the distribution of ET{(f)) are computed, 
by summing over the various cells {i} composing the PCAL, which is nearly symmetric in 
0: 

J2i Et cos(n0j) 

Vr, = ■■ 

J:^E^T 

Moments up to have been considered. Including the corresponding sine transforms, one 
can define two-dimensional vectors f„. The direction of vi for example, is an estimate of the 
direction of the reaction plane in a particular event. In principle then, the distributions of 
the relative phases of the Vn are meaningful too. 

A more transparent procedure might have been to fit an ellipse to the distribution Et{4>)- 
The center, the major and minor axes, and the orientation of the ellipse would provide a small 
set of physically meaningful parameters by which to characterise the distributions. The finite 
detector segmentation in limits the highest frequency f„ containing significant information. 
In addition, fluctuations are large on an event by event basis, and high frequency noise is 
present, further obscuring the meaning of higher moments. We consider it unlikely that much 
real information is present beyond vi and V2, and it is hard to imagine what physics could 
lead to such azimuthal distributions. In any case, only vi shows any significant systematic 
effects (dependence on centrality for example) in the present data. 

We have calculated vi and V2 from ARC events, simulating the segmentation of PCAL 
using sixteen bins in (p. The pseudorapidity interval considered is 1.85 < 1] < 4.5. 
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A comparison of the ARC moments with data is shown in Fig(|l]). One can see clearly 
that the distributions of f„ calculated from ARC agree quite well with the experimental ones, 
both in shape and in dependence on centrality. We plot only the magnitudes of f^, and the 
Jacobian factor (1/fn) is not included. The characters of the two dimensional distributions 
for Vi and V2 differ: the V2 distributions have a maximum at \v2\ =0 for all centralities. 
The distributions for vi are also peaked at zero for the most peripheral and most central 
event samples. But for semi-central events vi has a maximum away from zero. This is just 
the expected systematics of flow with centrality, the position of the maximum in vi being 
roughly related to the size of the transverse flow. 

Plots of the triple differential cross sections from ARC, for an event sample correspond- 
ing to the E877 central cut (10% of PCAL Et) are in Fig(^). For this cut, no flow is visible 
in the pion spectra. The pion distributions are symmetric around the beam axis, seemingly 
consistent with experiment 0. The proton spectra shift progressively towards positive {P^) 
with increasing rapidity. Additionally, at mid-rapidity the distribution is distorted, appear- 
ing roughly elliptical, with the major axis oriented along the {Py) direction. This anisotropy 
vanishes near beam rapidity, where one sees only the shift toward positive {Px)- This would 
appear to be the so-called 'squeeze out', in which protons are ejected symmetrically out of 
the reaction plane near mid-rapidity. 'Squeeze out' has not been not observed as yet at the 
AGS (it would be difficult at best to detect by analysing the but we note that quite 
a good definition of the reaction plane is probably necessary to discern it. Here, we are 
presenting essentially raw results from ARC. Folding these results with detector acceptances 
could possibly have a large effect. 

For the peripheral event sample Fig (^ the triple differential cross-sections are quite 
different. Most strikingly, the pion distributions show a clear shift toward negative (Px), 
i. e. opposite to the protons. This is the so-called 'anti-flow', and that it appears in the 
peripheral sample but is absent in the central one suggests that the effect is in large part 
due to shadowing by spectator nucleons. We note that E866 has detected an asymmetry in 
the pion distributions for semi-peripheral events [0]. 
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Our previous study of heavy ion data at BEVALAC energies ^ using ARC has shown 
that a good description of the observed flow can be obtained across a range of energies and 
projectiles. Central to achieving this agreement at low energy was a correct treatment of pp 
scattering, specifically the collisions leading to two-body final states. Kinetic pressure in the 
cascade was found to produce on the order of 90% of the observed fiow in the heaviest sys- 
tems, once angular momentum conservation was properly taken into account by preserving 
the reaction plane in such collisions. If the reaction plane were randomized instead, about 
70% of the fiow would be obtained. Only an additional 10% difference in the flow is ascribed 
to a potential-like effect in two-body scattering Proton rapidity and transverse mass 
spectra were also found to be in agreement with data and to be nearly independent of the 
choice of orbits and the randomisation of the reaction plane. The conclusion is that at these 
low energies little potential or mean fleld enters: the pure cascade accounts for the flow. 

Because the A^A^ interaction is mostly repulsive at the higher LBL energies, the classical 
cascade can represent scattering well using dominantly repulsive outgoing paths for two- 
body flnal states. In actual potential modeling the character of the classical orbits will 
vary with impact parameter and energy. One expects repulsive orbits for small impact 
parameters probing the inner core, and attractive orbits passing through the outer regions. 
At sufficiently high energy the description in terms of real potential scattering breaks down, 
and pion production becomes the dominant process. A reasonable guess is that somewhere 
above the threshhold for two pion production it is not sensible anymore to talk about 
potential models. In ARC, production of up to two pions pp nj^iT proceeds entirely through 
excitation of the AA intermediate state. Hence this scattering can still be described in terms 
of two baryon states. For flrst collisions at AGS energies the average pion number has risen 
to 3.8 per inelastic collision, the two body optical theorem is violated and no description in 
terms of potential scattering is possible. The potential, if any, is clearly dominated by its 
absorptive, imaginary part, and the scattering has become diffractive in nature. What then 
can one say about the nature of the classical orbits at such high energies? Diffraction from 
a black disk occurs in the absence of either repulsion or attraction in the potential, therefore 
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the best guess is probably that the orbits should be equally often repulsive or attractive. 

In constructing a cascade to be used at all energies between BEVALAC and AGS, one 
must interpolate between the mostly repulsive A^A^ scattering occuring at low energies and 
the diffractive A^A^ scattering taking place at higher energies. Heavy ion collisions at the 
AGS involve individual A^A^ scatterings taking place at all energies equal to or below the 
first collision energy. The simplest procedure which builds in the essential physics is to use 
repulsive orbits below some energy, Emm, a 50/50 mixture of repulsive and attractive orbits 
above a higher energy, Emax, and to interpolate smoothly between these limits. In practice 
we take Efjiin ^ 

300 MeV and E 

max ~ 500 MeV corresponding crudely to the effective 2,tt 
and Sir thresholds. We can distinguish this approach quite sharply from one involving mean 
fields: what we propose is not at all a medium dependent effect, rather it is an energy 
dependent treatment of the free space two body scattering. This is justified, we think, by 
the nature of the A^A^ interaction. 

In Fig(^) we display the effect on (P^.) of choosing an energy dependent mixture of 
repulsive and attractive orbits versus the 50/50 mixture. Choosing purely repulsive orbits 
at BEVALAC energies produced a rather small (< 10%) increase in flow discussed above. At 
AGS energies, however, the sensitivity to the classical orbits can be much larger, resulting in 
perhaps as much as 60% change in (P^) for a central cut. Indeed, the 'potential' one would 
need to turn all the orbits outward for p = 10.8 GeV/c would have to be very large indeed 
and completely unrealistic. If all the orbits are made attractive, interestingly, the flow is 
reduced to near zero. 

Summarising, present analysis of the very interesting E877 experiments on azimuthal 
asymmetry, combined with previous success at the BEVALAC, lends confidence to the cas- 
cade description of sideways flow. Certainly, the ARC calculations are in substantial agree- 
ment Fig.(0) with the higher energy data. The absolute flow magnitude is, however, more 
sensitive to the potential contribution at 11.6 GeV/c than at 2 GeV/c. The method of inter- 
polation between real potential and diffractive two body scattering is straightforward, but 
its full validation awaits the detailed E877 flow results and triple differential cross-sections. 
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Nevertheless, it is interesting to speculate on the predictions of this model at intermediate 
beam momenta, 2 — 8 GeV/c, relative to ongoing experiments with the EOS TPC Q], now 
installed at the AGS. At these, lower, AGS energies for Au + Au one may reasonably expect 
to form a largely equilibrated compound system, still at high baryon density. So, the likeli- 
hood of nucleating appreciable amounts of quark-gluon plasma should be optimized. The 
signature would be a considerable reduction in flow following from a reduction in kinetic 
pressure. The pioneering BEVALAC, E877 and E866 experiments provide needed guidance, 
fixing the flow at both low energy and at the full AGS Au energy, and setting the stage for 
perhaps dramatic results in between. 

The authors thank S.H.Kahana for useful discussions. This work has been supported by 
DOE grants No. DE-FG02-93 ER40768, DE-AC02-76 CH00016, and DE-FG02-92 ER40699. 
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FIG. 1. Moments of the distribution Et{4>). 
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FIG. 2. Proton and pion transverse momentum distributions relative to the event plane 
(x-direction) for 10% central events. Pion distributions are symmetric at all rapidities, 
indicating no pion anti-flow at this centrality. This seems to be consistent with the current 
experimental picture. Proton distributions show a clear shift towards positive {Px)- A 
spectator component is present in the proton distribution for the most forward rapidity. 
The deformation of the proton distributions along {Py) indicates proton squeezout occurs 
at mid-rapidity. 
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FIG. 3. Peripheral pion and proton momentum distributions (6 > 3 fm). Proton distri- 
butions appear symmetric at all rapidities: {Px) is small for peripheral samples. A large 
spectator component is present at the most forward rapidity. Pion distributions show a dis- 
tinct shift to negative {Px)'- that this pion 'anti-flow' appears only in the peripheral sample 
indicates it is largely due to shadowing. 
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FIG. 4. {Px) for protons and pfons. Central (10%) cut based on PCAL Et. Energy 
dependent mixture of repulsive scattering against 50/50 mixture at all energies. 
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This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9604008vl 
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